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We will describe a new superconductivity mechanism, proposed by the authors in [1], 
which is based on a topologically ordered ground state rather than on the usual Landau 
mechanism of spontaneous symmetry breaking. Contrary to anyon superconductivity it 
works in any dimension and it preserves P-and T-invariance. In particular we will discuss 
the low-energy effective field theory, what would be the Landau-Ginzburg formulation for 
conventional superconductors. 

1 Introduction 

For many years the theory of phase transitions was entirely understood in terms of the 
Landau-Ginzburg theory, based on symmetry breaking, order parameters and on the math- 
ematical framework of group theory. With the discovery of the fractional quantum Hall 
liquids [2], that are incompressible and exist only at some "magical" filling fractions, it was 
understood that the internal order characterizing these states is a new type of order, different 
from any other known type of order. This new type of order is the topological order [3], a 
particular type of quantum order 

Quantum order describes the zero-temperature properties of a quantum ground state, 
and characterizes universality classes of quantum states, described by complex ground state 
wave-functions [3]. Quantum phase transition are characterized by changes in the quan- 
tum entanglement properties of these complex ground state wave-functions. Topological 
order is a special type of quantum order whose hallmarks are the presence of a gap for all 
excitations (incompressibility) and the degeneracy of the ground state on manifolds with 
non-trivial topology [3]. In the case of the fractional quantum Hall effect, which is a (2-1-1)- 
dimensional system, another hallmark is the presence of excitations with fractional spin and 
statistics, called anyons [4]. The long-distance properties of these topological fluids can be 
explained by an infinite-dimensional VFi+oo dynamical symmetry [5], and are described by 
effective Chern-Simons field theories [6] with compact gauge group, which break P-and 
T-invariance. 

The topologically ordered superconductors we propose have a long-distance hydrody- 
namic action which can be entirely formulated in terms of generalized compact gauge fields. 
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the dominant term being the topological BF action [7]. The BF theory has many applica- 
tions, from 2D gravity [8] to the mathematical characterization of embedded manifolds [9]. 
It turns out that it is also the general model that describes the long-distance behaviour of 
systems that exhibit P- and T-preserving topological order in any dimension. It reduces 
to a doubled Chem-Simons model in (2-1-1 )-dimensions [10]. The existence of such non- 
conventional superconductors is also supported by purely algebraic considerations [11]. 

In [12], we have proven that planar Josephson junction arrays (JJA) can be exactly 
mapped onto an Abelian gauge theory with a mixed Chern-Simons term (BF-model): JJA 
provide thus a first concrete example of topological superconductors. The Abelian gauge 
theory exactly reproduces the phase diagram of JJA and the insulator/superconductor quan- 
tum phase transition at T = [13]. JJA have also been recently considered by several other 
authors [14, 15], as controllable devices which exhibit topological order. 

The idea that gauge fields can be used to model the long distance behaviour of con- 
densed matter systems was originally proposed in [17], and particularly exploited for pla- 
nar systems [18]. In a nutshell, the idea is that charge fluctuations around a given ground 
state are decribed by a conserved current j'^, which in (2-1-1) dimensions can be represented 
in terms of a gauge field according to j'* oc e^'^daBy. For a wide class of systems 
the effective action governing the dynamics of the charge fluctuations is quadratic in the 
gauge fields at long distances [17] . Clearly this effective action must also be gauge 
invariant, reflecting the original gauge invariance of the current j^: one obtains thus an 
effective gauge theory at long distances (which is not necessarily relativistic). The ground 
states of a wide class of planar condensed matter systems [19] can be classified according 
to the lowest derivative terms appearing in their effective gauge theories at long distances. 
This way, the Chern-Simons term describes incompressible quantum fluids (quantum Hall 
states) and chiral spin liquids [20] while the Maxwell term describes a (2-dim.) superfluid 
(superconductor) . 

The case of JJA is different in the sense that the gauge theory description is an exact 
mapping and not only an effective theory [12]: JJA can be exactly mapped onto an Abelian 
gauge theory with two gauge fields describing a current of charges and a current of vortices 
coupled by a mixed Chem-Simons term: 

C = -^,F,^F^^ + 1^^,6^-9.5. - ^/,./-^ , (1) 

Global superconductivity in planar JJA is thus the simplest example of the new mechanism 
of superconductivity we propose. 

We argue [21] that a less trivial example are frustrated JJA: these may support a topo- 
logically ordered superconducting ground state characterized by a non-trivial ground state 
degeneracy on the torus. These superconducting quantum fluids provide explicit examples 
of systems in which superconductivity arises from a topological mechanism rather than 
from the usual Landau-Ginzburg mechanism. 

The paper is organized as follows. In Section 2 we will review the definition of the 
Chern-Simons model on the lattice. In Section 3 we show how planar JJA can be ex- 
actly mapped onto a gauge theory with two Maxwell fields coupled by a mixed, periodic, 
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Chern-Simons term. Section 4 is devoted to the study of the phase diagram, showing how 
the superconductor/insulating transition of the JJA is reproduced by the gauge theory. In 
Section 5 we present the general theory of topological superconductors in any space-time 
dimension. Section 6 is, instead, devoted to the study of the role of frustration. 

2 Lattice Chern-Simons model 

Our model ([T]) can be rewritten in terms of the dual field strengths 

r Uu = d^B,-d,B^, (2) 

as follows (throughout this paper we use units such that c = 1 and h = 1.) 

For later convenience we have introduced a magnetic permeability rj, equal for the two 
gauge fields. The coupling constants and (7^ have dimension mass, whereas the co- 
efficient K of the mixed Chern-Simons term is dimensionless. Note that we take to 
represent a pseudovector gauge field, so that the mixed Chern-Simons term does not break 
the discrete symmetries of parity and time reversal. 

The action corresponding to (O is separately invariant under the two Abelian gauge 
transformations 

B^ -^B^ + d^Lo , (4) 

with gauge groups Ra and Rb, respectively. Moreover, the action is also invariant under 
the duality transformation 



A^ B^ , 



(5) 



so that the model is self-dual. 

The Lagrangian ^ can be easily diagonalized by the linear transformation 



Af, = J- (a^ + b^) , 

5/. =.[^{a^-b^). (6) 
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In terms of these new variables the model Q describes a free theory, 



eg \r] J 2tt eg \r] J 2tx 

(V) 

where G^ and g^ are the dual field strengths for the new gauge fields and 6^, respectively. 
This Lagrangian describes a doublet of excitations with topological mass [22] 



\K,\eq 

m = 4^, (8) 



and spectrum 



£;(q) = Jm2 + i|q|2 . (9) 



In the following we shall formulate a Euclidean lattice version of the above Chem- 
Simons model. To this end we introduce a three-dimensional rectangular lattice with lattice 
spacings in the three directions. In particular we shall take the lattice spacings li = I2 = I 
and identify Iq with the spacing in the EucUdean time direction. Lattice sites are denoted 
by the three-dimensional vector x; the gauge fields A^{x) and B^{x) are associated with 
the links (x, /i) between the sites x and x + fi, where fi denotes a unit vector in direction /i 
on the lattice. 

On the lattice we introduce the following forward and backward derivatives and shift 
operators: 

d,,f{x) = ^ , S^f{x) = fix + l^n) , 

d,f{x) ^M^Ii^^Ml^ s,f{x)^f{x-l,f^). (10) 

Summation by parts on the lattice interchanges both the two derivatives (with a minus sign) 
and the two shift operators; gauge transformations are defined using the forward lattice 
derivative. Corresponding to the two derivatives in [TO] , we can define also two lattice 
analogues of the Chem-Simons operators e^auda [23] [12] : 

where no summation is implied over equal indices ji and v. Summation by parts on the 
lattice interchanges also these two operators (without an extra minus sign). The operators 
[TTl are both local and gauge invariant, in the sense that 

k^udu — d^k^i/ = , k^ydy = d^k^i^ = , (12) 

and their product reproduces the relativistic, EucUdean lattice Maxwell operator: 

kfxakau — kfj^akau — ^fj,u^ ~l~ d^dy , (13) 
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where = (i/^d^ is the three-dimensional Laplace operator. Using kf^^ we can also define 
the lattice dual field strengths as 

= kf^uB, . (14) 

The identity (fT3] ) then tells us that we can simply write the relativistic, Euclidean lattice 

Maxwell terms as J2x -^m-^/^ f^lf^l■ 

Using all these definitions we can now write the Euclidean lattice partition function of 
our model Q as follows: 

Zcs = |^^^^| eM-Scs) , 

Scs = E ^ (^^oFo + F.i^.) - ^'-^A,k,.B^ + |! (i/o/o + (15) 

where we have introduced the notation VA^ = Y\[x,fj.) dA^{x) and gauge fixing is under- 
stood. 

For later convenience we introduce also the finite difference operators 

= l^dfj, , = l^dfj. , (16) 

where no summation over equal indices is implied. Correspondingly, we introduce also the 
finite difference analogue of the operators k^iy and /c^i^: 

= S^S^au^a 5 = d^Q^yAoiSu . (T7) 

These satisfy equations analogous to (fT2l) and ([T3] ) with all derivatives substituted by finite 
differences. 



3 Josephson junction arrays 

JJA we analize are quadratic, planar arrays of spacing I of superconducting islands with 
nearest neighbours Josephson couplings of strength Ej. Each island has a capacitance Cq 
to the ground; moreover there are also nearest neighbours capacitances C. The Hamiltonian 
characterizing such systems is thus given by 

^ = E + E (? (^y - ^x)' + Ej{l- cos N (a>y - CD J)) , (18) 

X ^ <xy> \ ^ / 

where boldface characters denote the sites of the two-dimensional array, < xy > indicates 
nearest neighbours, 14; is the electric potential of the island at x and <I>x the phase of its 
order parameter. For generality we allow for any integer N in the Josephson coupling, so 
that the phase has periodicity I-k/N: obviously = 2 for the real systems. 
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With the notation introduced in the previous section the Hamiltonian (fTSl ) can be rewrit- 
ten as 

H = J2lyiCo-CA)V + J2 Ej (1 - COS N (A,$)) , (19) 

X x,j 

where A = AjAi is the two-dimensional finite difference Laplacian and we have omitted 
the explicit location indeces on the variables V and <I>. 

The phases <l>x are quantum-mechanically conjugated to the charges Qx on the islands: 
these are quantized in integer multiples of N (Cooper pairs for N = 2): 

Q = qeNpo , poeZ (20) 

where is the electron charge. The Hamiltonian ( fT9l ) can be expressed in terms of charges 
and phases by noting that the electric potentials are determined by the charges Qx via a 
discrete version of Poisson's equation: 

[Co - CA) ^x = Qx . (21) 

Using this in ( [T9] ) we get 

H = J2 N^Ec POg^^Po + J2 Ejil- cos N (A,«>)) , (22) 

X C x,i 

where Ec = q^/2C. The integer charges po interact via a two-dimensional Yukawa poten- 
tial of mass ^JCqJC jl. In the nearest-neighbours capacitance limit C » Co, which is ac- 
cessible experimentally, this becomes essentially a two-dimensional Coulomb law. ^From 
now on we shall consider the limiting case Co = 0. In this case the charging energy Eq 
and the Josephson coupling Ej are the two relevant energy scales in the problem. These 
two massive parameters can also be traded for one massive parameter \/2N'^EcEj, which 
represents the Josephson plasma frequency and one massless parameter Ej/Ec- 

The zero-temperature partition function of the Josephson junction array admits a 
(phase-space) path-integral representation [24]. Since the variables po are integers, the 
imaginary-time integration has to be performed stepwise; we introduce therefore a lattice 
spacing Iq also in the imaginary-time direction. This has to be just smaller of the typical 
time scale on which the integers po vary, in the present case the inverse of the Josephson 
plasma frequency: Iq < O {l/^/2N'^EcEj\ . We thus get the following partition function: 



-PO + E ^o^J (1 - cos N (A,;$)) , (23) 

where now the sum in the action S extends over the three-dimensional lattice with spacing 
Iq in the imaginary time direction and I in the spatial directions. 



j-+iv/N 

S =J2 PoAo^ + N'^Ech Po- 
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In the next step we introduce vortex degrees of freedom by replacing the Josephson 
term by its Villain form: 



exp(-5) , 

{"J 



(24) 



Strictly speaking, this substitution is valid only for IqEj » 1; however the Villain approx- 
imation retains all most relevant features of the Josephson coupling for the whole range of 
values of the coupling Ej and therefore we shall henceforth adopt it. 

We now represent the Villain term as a Gaussian integral over real variables pi and 
we transform also po to a real variable by introducing new integers vq via the Poisson 
summation formula 

fc=+oo n=+oo 

exp{i2TTkz) = 6{z-n) . (25) 

k=—oo n=-oo 

By grouping together the real and integer p and v variables into three-vectors and v^, 
^ = 0, 1, 2 we can write the partition function as 



r p+tt/N 

Z = E / exp(-5) , 

,1 J J -TV IN 



/ 27r \ 1 

S =Y -iNp^ (^A^cD + —v^j + N^Eck Po^Po + • (26) 

Following [25] we use the longitudinal part of the integer vector field to shift the 
integration domain of To this end we decompose Vfj, as follows: 

= A^m + A^a + K^^ip^ , (27) 

where m G Z, |a| < 1 and K^i^ defined in ([17] ) . Here the vectors tp^ are not integer, but 
they are nonetheless restricted by the fact that the combinations = K^^Vy = K^^Kau'^u 
must be integers. The original sum over the three independent integers {w^} can thus 
be traded for a sum over the four integers {m, q^} subject to the constraint A^g^ = 0. 
The sum over the integers {m} can then be used to shift the $ integration domain from 
[—tt/N^+tt/N) to (— oo,+oo). The integration over <I> is now trivial and enforces the 
constraint A^p^ = 0: 



S =Y -^2^ Pf^Kf^ui^u + N\Ec Po^Po + ■ (28) 
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We now solve the two constraints by introducing a real gauge field and an integer 
gauge field o^: 

Pfi = K^yhu , bfj,£ R, 

Qfi = K^iyttu , a^e Z . (29) 



By inserting the first of these two equations and by summing by parts, the first term in 
the action ( [281 ) reduces to J2x "^^vrft^g^. By inserting the second of the above equations 
and by summing by parts again, this term of the action finally reduces to the mixed Chern- 
Simons coupling J2x — ^27r a^K^yhy. Using the Poisson formula (|25] ) we can finally make 
also real at the expense of introducing a set of integer link variables {Q^} satisfying the 
constraint A^Q^, which guarantees gauge invariance: 



Wm} 

1 iD^ 

S =y2-i2Tia^K^,yK + N\EcPi)—rPo + ^;f^ + i2Tia^Q^. (30) 
^ -A MqEj 



In this representation K^^hy represents the conserved three-current of charges, while 
K^ytty represents the conserved three-current of vortices. Note that, actually, both these 
conserved currents are integers (the factors of are explicit): indeed, the summation over 
{Qfi} makes a^j (and therefore also K^yUy) an integer, and then the summation over {a^} 
makes K^yhy an integer. The third term in the action (l30l ) contains two parts: the longi- 
tudinal part [pi^ describes the Josephson currents and represents a kinetic term for the 

charges; the transverse part {pj^ can be rewritten as a Coulomb interaction term for the 
vortex density go by solving the Gauss law enforced by the Lagrange multiplier feo- 

The partition function (l30l ) displays a high degree of symmetry between the charge and 
the vortex degrees of freedom. The only term which breaks this symmetry (apart from the 
integers Q^^) is encoded in the kinetic term for the charges (Josephson currents). This near- 
duality between charges and vortices has already been often invoked in the literature [16] 
to explain the experimental quantum phase diagram at very low temperatures. Here we 
introduce what we call the self-dual approximation of Josephson junction arrays. This 
consists in adding to the action in ([30l ) a bare kinetic term for the vortices (note that such 
a kinetic term is anyhow induced by integrating out the charge de|rees of freedom) and 
combining this with the Coulomb term for the charges into Yl,x n'%Ec "^^^ coefficient 
is chosen so that the transverse part of this term reproduces exactly the Coulomb term 
for the charges upon solving the Gauss law enforced by the Lagrange multiplier oq. The 
longitudinal part, instead, represents the additional bare kinetic term for the vortices. Given 
that now the gauge field has acquired a kinetic term, we are also forced to introduce new 
integers via the Poisson formula to guarantee that the charge current K^j^yby remains an 
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integer: 

ZsD = J2 J ^'^A' / ^^A* exp(-55Z)) , 



{Q^} 



Pi , 7r2g2 



SsD = Y] -i27r a^i^^^6^ + -p— + + i27ra^Q^ + i27rb^M^ , (31) 

where the new integers satisfy the constraint A^M^ = to guarantee gauge invariance. 
After a rescaling 

^0 = f-T, 0,0 , Ai = r— ai , 



_ 27r _ 27r 

-Bo = r— bo , Bi = —7=-bi , (32) 

VNlo VNl 



we obtain finally 

ZsD = E / '^\J -DB^expi-SsD) , 



(^oQo^o + IQiAi) + {IoMqBo + IMiBi) , (33) 
where Fi and fi are defined in (fT4l) and 

This is exactly the partition function of our lattice Chern-Simons model ( fT5] ) in the limit 
of infinite magnetic permeability r] = oo and with additional, integer-valued link variables 
and coupled to the two gauge fields. Note that, with the above identifications, the 
topological Chern-Simons mass ^ coincides with the Josephson plasma frequency: 



m 



V2N^EcEj . (35) 



In the physical case N = 2 this reduces to m = y/SEcEj. From the kinetic terms in ([33 
we can also read off the charge and vortex masses: 





1 


N 


rriq 




iirH^Ej 




1 


1 




Z2e2 


2NPEc 



(36) 

In the regime ml < 0{1), which is typically experimentally relevant, we can choose lo = I- 
in this case the infinite magnetic permeability constitutes the only non-relativistic effect in 
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the physics of Josephson junction arrays in the self-dual approximation. However, we ex- 
pect this non-relativistic effect to be irrelevant as far as the phase structure and the charge- 
vorticity assignements are concerned. Therefore, for simplicity, we shall henceforth con- 
sider the relativistic model, by setting Iq = I and rj = 1, although it is not hard to incorporate 
a generic value of r] into our subsequent formalism: 



ZsD = J2 J^\J^B^^M-Ssd) 



{Mm} 



/3 ^3^ ^3 



X 



Josephson junction arrays in the self-dual approximation constitute thus a further, exper- 
imentally accessible example of the ideas presented in [17] and [18]. The action in ( [37] ) 
provides in fact a pure gauge theory representation of a model of interacting charges and 
vortices, represented by the conserved currents 

3 

charge = ^ U R 



1 vortex 



= r k^yAy , (38) 

27rK2 



where the prefactors are chosen so that the quantum of charge is given by k, while the 
quantum of vorticity is given by 1/k (factors of and 2tt are absorbed in the definitions of 
the gauge fields and the coupling constants). 

In this framework, the mixed Chem-Simons term represents both the Lorentz force 
caused by vortices on charges (coupling of q'^^'^s^ "electric" gauge field A^) and, by 
a summation by parts, the Magnus force caused by charges on vortices (coupling of 0™'^*^^ 
to the "magnetic" gauge field B^). The integer-valued link variables and represent 
the (Euclidean) topological excitations [25] in the model. They satisfy the constraints 

d^,Qfl = , 

= . (39) 

In a dilute phase they constitute closed electric (Q^) and magnetic (M^) loops on the lattice; 
in a dense phase there is the additional possibility of infinitely long strings. Due to the 
constraints (l39l ) we can choose to represent these topological excitations as 

=lk^yXy, X^eZ, (40) 

and reabsorb them in the mixed Chem-Simons term as follows: 

Sso = ^...-^{a,- V {b. - f^y,) + . . . . (4.) 
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In this representation it is clear that the topological excitations render the charge-vortex 
coupling periodic under the shifts 



27r 

-^B^ + j-^b^, b^£Z. (42) 



In physical terms, the topological excitations implement the well-known [19] periodicity 
of the charge dynamics under the addition of an integer multiple of the flux quantum 1/k 
per plaquette and the (less-known) periodicity of the vortex dynamics under the addition of 
an integer multiple of the charge quantum k per site. 

If we would require that the full action (including charge-charge and vortex-vortex in- 
teractions) (l37l) be periodic under the shifts (l42l ) , then we would obtain the compact Chem- 
Simons model studied in [12] . In this case the relevant topological excitations would be 
essentially iX^ and iY/^^: since these can also describe finite open strings, there is the ad- 
ditional possibility of electric and magnetic monopoles [25] . As we showed in [12], these 
monopoles play a crucial role in the regime ml <^ I. 



4 Phase structure analysis 

In this section we investigate symmetry aspects and non-perturbative features of the model 
(|37] ) due to the periodicity of the charge-vortex interactions encoded in the mixed Chern- 
Simons term. As expected, these depend entirely on the topological excitations which en- 
force the periodicity. 

Upon a Gaussian integration the partition function (l37l) factorizes readily as 



ZsD = Zcs ■ Ztop , (43) 
where Zcs is the pure gauge part defined in ([T5] ) and 



Ztop = exp(-5Top) , 

{Qm> 

with m defined in dS]), describes the contribution due to the topological excitations. The 
phase structure of our model is thus determined by the statistical mechanics of a coupled 
gas of closed or infinitely long electric and magnetic strings with short-range Yukawa in- 
teractions. The scale represents the width of these strings. In our case it is of the 
order of the lattice spacing /. The third term in the action (|44l) . describing the topological 
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Aharonov-Bohm interaction of electric and magnetic strings, vanishes for strings separated 
by distances much bigger than (1/m): in this case the denominator reduces to m^V^ and, 
by using either one of the two equations in (l40l ) and the constraints ( [39] ) one recognizes 
immediately that the whole term in the action reduces to (ilir integer), which is equivalent 
to (this reflects the fact that the original charges and vortices satisfy the Dirac quantization 
condition). 

4.1 Free energy arguments 

In order to establish the phase diagram of our model we use the free energy arguments for 
strings introduced in [26]. 

The usual argument for strings with Coulomb interactions [26] is that interactions be- 
tween strings are unimportant for the phase structure because small strings interact via 
short-range dipole interactions, while large strings have most of their multipole moments 
canceled by fluctuations. This argument is even stronger in our case, where the interac- 
tion is anyway short-range. Therefore one retains only the self-energy of strings, which is 
proportional to their length, and phase transitions from dilute to dense phases appear when 
the entropy of large strings, also proportional to their length, overwhelms the self-energy. 
We shall also neglect the interaction term between electric and magnetic strings (imaginary 
term in the action ([44|) ). This is clearly a good approximation if both types of topological 
excitations are dilute. 

Thus, one assigns a free energy 

F = (^^G{ml) Q2 + l^G{ml)M^ - N (45) 

to a string of length L = IN carrying electric and magnetic quantum numbers Q and M, 
respectively. Here G{ml) is the diagonal element of the lattice kernel G{x — y) representing 
the inverse of the operator [m? — V^). Clearly G{ml) is a function of the dimensionless 
parameter ml. The last term in ( [45]) represents the entropy of the string: the parameter ^ 
is given roughly by /x = ln5, since at each step the string can choose between 5 different 
directions. In ( [45] ) we have neglected all subdominant functions of N, like a \n.N correction 
to the entropy. 

The condition for condensation of topological excitations is obtained by minimizing the 
free energy ([45] ) as a function of N. If the coefficient of in ([45] ) is positive, the minimum 
is obtained for A^ = and topological excitations are suppressed. If, instead, the same 
coefficient is negative, the minimum is obtained for A^ = oo and the system will favour 
the formation of large closed loops and infinitely long strings. Topological excitations with 
quantum numbers Q and M condense therefore if 

le^nGiml) Ig^KGiml) ^ 

^ — - + — ^ — - < 1 . (46) 

2/i 2n 

If two or more condensations are allowed by this condition one has to choose the one with 
the lowest free energy. 
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The condition (l46l ) describes the interior of an ellipse with semi-axes 2ii / {le^ KG{ml)) 
and 2fj,/{lg'^KG{ml)) on a square lattice of integer electric and magnetic charges. The 
phase diagram is obtained by investigating which points of the integer lattice lie inside the 
ellipse as its semi-axes are varied. We find it convenient to present the results in terms of 
the dimensionless parameters Im and e/g: 

mlG{ml)7r i g ^ ^ ' electric condensation , 

^ ^ ~^ 1 I > 1 , magnetic condensation , 

electric condensation , 




. 3 mlG{ml)'rr ' 

mlG(ml)TT ^ a ^ e ^ mlG(ml)n , .■ , .^^ 

— — — > 1 ^ < ,J7 n < - < 5^ — — , no condensation , (47) 

^ ^ I mlG(ml)Tv g /i ' ' ' / 

I e ^ miG(rni)n ^ magnetic condensation . 

As expected, these condensation patterns are symmetric around the the point e/g = 1, 
reflecting the self-duality of the model. In first approximation the electric (magnetic) con- 
densation phase is characterized by the fact that {Qfi} ({A/^}) fluctuate freely, while all 
Mf_i = (Q^ = 0). Within this approximation it is clearly consistent to neglect altogether 
the interaction term between electric and magnetic strings in (l45l) . Taking into account small 
loop corrections in the various phases can lead to a renormalization of coupling constants 
and masses and, correspondingly, to a shift of the critical couplings (m/)crit and (e/(7)crit 
for the phase transitions. A notable exception is the case in which there is only one phase 
transition: in this case the critical coupling is {e/g)cvit = 1 due to self-duality. 



4.2 Wilson and 't Hooft loops 

In order to distinguish the various phases we introduce the typical order parameters of lattice 
gauge theories [25], namely the Wilson loop for an electric charge q and the 't Hooft loop 
for a vortex 0: 




where and vanish everywhere but on the links of the closed loops, where they take 
the value 1 . Since the loops are closed they satisfy 

df,q^ = d^(l)t, = . (49) 

The expectation values (Lw) and {Lh) can be used to characterize the various phases. 
First of all they measure the interaction potential between static, external test charges q and 
—q and vortices (j) and —cj), respectively [25] . Secondly, by representing the closed loops 
q^ and 0^ as 

Qfj, = I'kfMu-^'i/ ) 

cP^ = Ik^^At , (50) 
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we can rewrite the Wilson and 't Hooft loops as 

Lh =exp|^mi 0X^/2^^/^^ , (51) 

which is a lattice version of Stoke's theorem, the integers and A'^ (= 0, ±1) representing 
the area elements of the surfaces spanned by the closed loops. The second terms of the 
expansions of {Lw) and {Lh) in powers of q and (j) measure therefore the gauge invariant 
correlation functions {F^{x)F^{y)) and {f^{x)f^{y)). Third, if we represent (f)^ as 



/2 ^ 

— k^ijAjj ■ , (52) 



we can also rewrite the 't Hooft loop as 



With the interpretation of as an external electromagnetic gauge potential the expecta- 
tion value of the 't Hooft loop measures the electromagnetic response of the system in the 
various phases. An analogous relation clearly holds for the Wilson loop. 

The expectation values of the Wilson and 't Hooft loops are easily obtained by combin- 
ing the definitions |48] withl37l 



{Lw) 



^"^^^ = Ztop(Q.,M,) ' ^''^ 

where the notation is self-explanatory. In the following we shall analyze these expressions 
in the various phases obtained in (|47] ). We shall mostly only indicate the form of small loop 
corrections: a full renormalization group analysis is beyond the scope of the present paper 
and we won't be able to predict the orders of the phase transitions. 

Let us begin with the electric condensation phase. In this phase the ground state con- 
tains many infinitely long electric strings Q^. These have a crucial effect on the gauge 
symmetry associated with the gauge field A^. To see this let us consider a gauge transfor- 
mation Afj^ — > Afj^ + d^A, where, for simplicity, we take A as a function of the component 

only. If we choose the usual boundary conditions = = at infinity, the change of 
the action (|37] ) under the above gauge transformation is given by 



^SsD = (a(x^ = +oo)Qi{x^ = +oo) - A(x^ = -oo)Qi{x^ = -oo 



(55) 
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In a dilute phase, with only small closed loops, Qi{x^ = +00) = Qi{x^ = —00) = 
and the action is automatically gauge invariant. In a dense phase, with many infinitely long 
strings, Qi{x^ = +00) and Qi{x^ = —00) are generically different from zero. Gauge 
invariance requires that ASsd vanishes modulo z27r. In the dense phase this is realized 
only if A takes the values 

A = —=n , ne Z , (56) 



at infinity. This means that, in the electric condensation phase, the global gauge symmetry 
is spontaneously broken down to the discrete gauge group Z, so that the total (global) 
symmetry of this phase is x Rq. 

The Wilson loop expectation value takes a particularly simple form if the external test 
charges are multiples of the charge quantum: q = nn, n G Z. In fact, since we sum over 
{Qjj}, the integer loop variables nq^ can be absorbed by a redefinition of the appropriate 
Q^'s, with the result 

{Lw{q = nK)) = 1 . (57) 

This indicates that, in this phase, external test charges q = nn are perfectly screened by 
the topological excitations and behave thus freely. In order to compute the Wilson loop 
expectation value for generic q we have to perform explicitly the sum over {Q^}. To this 
end we have to remember the constraint d^^Q^ = 0. We solve this constraint by representing 
QpL = Ikf^yUy and summing over {n^}, with the appropriate gauge fixing. We then use 
Poisson's formula [25] to turn this sum into an integral, by introducing a new set of integer 
link variables {k^} satisfying d^k^ = in order to guarantee the gauge invariance under 
iT-fj. ^ + Id^i. At this point we can perform explicitly the Gaussian integration over 
{n^}. In the approximation of neglecting terms proportional to V^/m^ (keeping such 
terms would not alter substantially the result) the new integers {A:^} can be absorbed by a 
redefinition of the magnetic topological excitations {M^}, giving the result: 

-^corr. 



w 



Zcovv. {Qim — 0) 

21 ^-V2" 



^corr.M = E exp^(-^M^^M, + i27r^yl^Mj . (58) 



{Mpjloops X 

Since the magnetic topological excitations are in a dilute phase we have to sum only over 
small closed loops: in this phase the dominant part of ln{L\v) vanishes for generic q and 
the whole result is given by small loop corrections. These are identical in form to the 
small loop corrections for the correlation functions in the low-temperature phase of the 
three-dimensional XY model; correspondingly the Wilson loop expectation value can be 
computed by exactly the same low-temperature expansion used for the XY model. The first- 
order term in this expansion is obtained by considering only the smallest possible lattice 
loops and gives the result 



kI 



{Lw) = exp I 2e " e ^ 



COS ( 27r-g^ ) - 1 



(59) 
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The periodicity of this result is a direct consequence of the spontaneous symmetry breaking 
Ra —>■ Za- This implies also that the gauge invariant correlation function reduces to 

{F^ix)F,{y)) oc (5^,V2 - d^d,) ^ , (60) 

which is essentially a contact term on the scale of the lattice spacing. 

The computation of the 't Hooft loop expectation value follows exactly the same lines 
as the above computation of the Wilson loop. The results is 

[Lh) = exp -— ^ ^ ^ 1 



Zcorr (0^) = E exp ( M^^M, + 2k<P M^^c/)^ ] . (61) 

{A/;j}loops 



X 



The first few terms in the expansion of the small loop corrections can again be computed 
with the same techniques as in the low-temperature phase of the XY model. One finds that 
their contribution amounts to perturbative corrections of the Coulomb coupling constant 
g^fC'cfP' /2l of the dominant term in (|6T]) . 

^From (|6TI ) we can extract the nature of the electric condensation phase. First of all, by 
considering, as usual, a rectangular loop of length T in the imaginary time direction and 
of length R in one of the spatial directions and computing the dominant large-T behaviour 
of ln(L//) we find that the interaction potential between external test vortices of strength 
and — (/) is proportional to \nR. Vortices are thus logarithmically confined, which amounts 
to the Meissner effect. Secondly, by using the representations (l50l ) and |5T] , we find the 
correlation function 

{M-)fAy)) oc '^^^^ ^ , (62) 

which is long-range, indicating that the "S^-photon" is massless. This is the massless ex- 
citation associated with the spontaneous symmetry breaking of the global gauge symmetry 
Ra —>■ Za- Third, by using the representations (l52l) and (l53l) . we find that the induced 
electromagnetic current is given by 

jr « - ^) ^r- , (63) 

which is the standard London form. We thus conclude that the electric condensation phase 
is actually a superconducting phase. 

No further computation is needed to extract the nature of the magnetic condensation 
phase: this is the exact dual of the electric condensation phase just described. Specifically, 
the global gauge symmetry associated with is spontaneoulsy broken down to Zb, so 
that the total symmetry of this phase is Ra x Zb- Correspondingly, the "^^-photon" is 
massless and the {F^{x)Fi,{y)) correlation function is long-range. Electric charges are log- 
arithmically confined, which means that an infinite energy (voltage) is required to separate 
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a neutral pair of charges. We call this phase with infinite resistance a superinsulator. In real 
Josephson junction arrays we expect however the conduction gap to be large but finite due 
to the small ground capacity Co, resulting in a normal insulator. 

If mlG{ml)ir/ fi > 1 a third phase can open up between the superconducting and su- 
perinsulating phases. In this third phase both the electric and the magnetic topological 
excitations are dilute. Far away from the phase transitions and to first approximation we 
can neglect them altogether. This gives the result 



{Lw) = exp -— — 2 ^"^72 

\ 21k — 



{Lh) = exp y-^^f- 0M ^2'^rv2 '^'^J • ^^"^^ 

Small loop corrections to these results can be obtained by restricting the {Q^} and {M^} 
sums in (l54l ) to small closed loops and using again the same techniques as in the low- 
temperature expansion of the XY model. These will lead to perturbative corrections of 
the coupling constants and masses in (l64l ): however the first-order result (l64l ) is enough to 
establish the nature of this phase. The global symmetry characterizing this phase is Ra x 
Rb and, corrrespondingly, both "photons" are massive, resulting in short-range correlation 
functions {F^{x)Fy{y)) and {f^{x)f,y{y)). Both charges and vortices interact via short- 
range Yukawa potentials and behave thus freely when separated by distances larger then the 
scale (1/m). In presence of any dissipation mechanism (which would not alter the other 
two phases) this third phase corresponds thus to a metallic phase of the Josephson junction 
array. 

In conclusion we can represent the phase diagram of our model as follows: 
mlG{ml)iT ^ ^ I I < 1 , superconductor {Z^Rb) , 



H I I > 1 , superinsulator {Ra x Zb) , 



r f < ^iJl^i)^ ^ superconductor {Za x Rb) , 

mlGimljir „ ^ e ^ mlG(ml)Tr ^ i / n r> \ 

y-^ > 1 micU)- < f < — TT^ ' "'^^^^ ' ^65) 

1^ I ^ miGirni)!^ ^ supcriusulator {Ra x Zb) , 

where we have indicated in parenthesis the global symmetries of the various phases. 

A numerical computation of the function mlG{'ml)'ir/ ^ for the value fi = ln5, gives 
an indication that a window for the metallic phase is open for ml just larger than 1 , while 
in the regime ml < 0(1), relevant for Josephson junction arrays, a single phase transition 
from a superconductor to a superinsulator at {e/g) = I is favoured. 

The experimental results for Josephson junction arrays, which are essentially resistance 
measurements as a function of temperature in arrays with 0(10'^) cells, shows, etrapolating 
at zero-temperature, a quantum phase transition between an insulator and a superconductor 
in the vicinity of the self-dual point Ej/Ec = ~ 0.2. This is in perfect agreement 
with (e/g) = 1. In fact, we have: Ej/Ec = (2g'^) / {e^ix'^) and, for {e/g) = 1, we get 
Ej/Ec = 2/7r2 ~ 0.2. 
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5 Superconductor with topological order 

Topologically ordered superconductors have a long-distance hydrodynamic action which 
can be entirely formulated in terms of generaUzed compact gauge fields, the dominant term 
being the topological BF action. 

BF theories are topological theories that can be defined on manifolds M^+i of any 
dimension (here d is the number of spatial dimensions) and play a crucial role in models of 
two-dimensional gravity [8]. In [12] we have shown that the BF term also plays a crucial 
role in the physics of Josephson junction arrays. 

The BF term [7] is the wedge product of a p-form B and the curvature dA of a (d-p) 
form A: 

k f 

where A; is a dimensionless coupling constant. This can also be written as 

Sbf = ^ [ Ad-p A dBp . (67) 
27r JMd+i 

The integration by parts does not imply any surface term since we will concentrate on 
compact spatial manifolds without boundaries and we require that the fields go to pure 
gauge configuration at infinity in the time direction. Indeed this action has a generahzed 
Abehan gauge symmetry under the transformation 

B^B + ri, (68) 

where 77 is a closed p form: drj = 0. Gauge transformations: 



A^A + ^, (69) 

with ^ a closed (d-p) form instead, change the action by a surface term. This, however 
vanishes with the boundary conditions we have chosen. 

Here we will be interested in the special case where Ai is a 1-form and, correspond- 
ingly, Bd-i is a (d-l)-form: 

Sbf = ^ [ AiA dBd-i . (70) 

27r JMd+i 

In the special case of (3+1) dimensions, B is the well-known Kalb-Ramond tensor field 
B^, [27]. 

In the application to superconductivity, the conserved current ji = *dBd-\ represents 
the charge fluctuations, while the generalized current jd^i = *dAi describes the conserved 
fluctuations of (d-2)-dimensional vortex lines. As a consequence, the form B^-i must be 
considered as a pseudo-tensor, while Ai is a vector, as usual. The BF coupling is thus P- 
and T-invariant. 
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The low-energy effective theory of the superconductor can be entirely expressed in 
terms of the generalized gauge fields Ai and B^-i- The dominant term at long distances is 
the BF term; the next terms in the derivative expansion of the effective theory are the kinetic 
terms for the two gauge fields (for simplicity of presentation we shall assume relativistic 
invariance), giving: 



Stm = Jm, ., ^dAi A *dAi + ^Ai A dBd 



-1 



+ t^dBd-i A *dBd-i , (71) 

where and g"^ are coupling constants of dimension m~'^+^ and m'^^^ respectively. 

The BF-term is the generalization to any number of dimensions of the Chern-Simons 
mechanism for the topological mass [22]. To see this let us now compute the equation of 
motion for the two forms A and B: 



-ird * dBd-i = — 
5 27r 

and 



2 d * dBd_i = —dAi , (72) 



-^d *dAi = — 
Applying d* on both sides of (1721) and (1731) we obtain 



-2 d * dAi = —dBd-i ■ (73) 



d*d* dAi - ^d * dBd^i = , 

d*d*dBd-i- ^d*dAi=0 . (74) 

The expression *d* is proportional to 6, the adjoint of the exterior derivative [28]. Substi- 
tuing d * dBfi^i and d * dAi in ( 1741 ) with the expression coming from ( 1721 ) and ( 1731 ) we 
obtain 



(A + m^) dAi = , 

(A + m2) dBd-i = , (75) 

where A = ci5 (when acting on an exact form) and m = is the topological mass. 
This topological mass plays the role of the gap characterizing the superconducting ground 
state. Note that the gap arises here from a topological mechanism and not from a local 
order parameter acquiring a vacuum expectation value. Equations (1721 ) and (1731 ) tell us 
that charges are sources for vortex line currents encircling them and viceversa. This is the 
coupling between charges and vortices at the origin of the gap. A related mechanism for 
topological mass generation in (3+l)-dimensional gauge theories is the generalization of 
the Schwinger mechanism proposed in [29]. 
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Let us now consider the special case of (2+1) dimensions (d=2). In this case also B 
becomes a (pseudo-vector) 1-form and, correspondingly the BF term reduces to a mixed 
Chem-Simons term. This can be diagonalized by a transformation A = B = a — b, 
giving 



The result is a doubled Chem-Simons model for gauge fields of opposite chirality. This 
action, including its non-Abelian generalization with kinetic terms was first considered 
in [30]. It is the simplest example of the class of P- and T-invariant topological phases 
of strongly correlated (2-i-l)-dimensional electron systems considered in [10]. Indeed, the 
BF term is the natural generalization of such doubled Chem-Simons models to any di- 
mension. Doubled (or mixed) Chem-Simons models are thus particular examples in two 
spatial dimensions of a wider class of P- and T-invariant topological fluids that have a su- 
perconducting phase. These fluids are described by the topological BF theory with compact 
support for both gauge fields. 

Topological BF models provide also a generalization of anyons to arbitrary dimensions. 
While in (2+1) dimensions fractional statistics arises from the representations of the braid 
group, encoding the exchange of particles, in (3+1) dimensions it arises from the adiabatic 
transport of particles around vortex strings and, in (d+1) dimensions, from the motion of 
an hypersurface S/j around another hypersurface S^-Zi. The relevant group in this case is 
the motion group and the statistical parameter is given by ^h{d — h), where k is the BF 
coupling constant [31]. 

Let us now illustrate the mechanism of superconductivity. To this end we shall from 
now on consider only rational k = ^ with ki integers, and specialize to manifolds M^+i = 
Md X whith Ri representing the time direction. 

The compactness of the gauge fields allows for the presence of topological defects, both 
electric and magnetic. The electric topological defects couple to the form Ai and are string- 
like objects described by a singular closed 1-form Qi. Magnetic topological defects couple 
to the form B^-i and are closed (d-l)-branes described by a singular (d-l)-dimensional 
form Qd-i- These forms represent the singular parts of the field strenghts dAi and dB^-i, 
allowed by the compactness of the gauge symmetries [25], and are such that the integral 
of their Hodge dual over any hypersurface of dimensions d and 2, respectively, is 27r times 
an integer as can be easily derived using a lattice regularization. Contrary to the currents 
j\ and jd-i, which represent charge- and vortex-density waves, the topological defects 
describe localized charges and vortices. In the effective theory these have stmcture on the 
scale of the ultraviolet cutoff. 

We will not discuss here the conditions for the condensation of topological defects, but 
we will show, instead that the phase of electric condensation describes a superconducting 
phase in any dimension. A detailed analysis would require the use of an ultraviolet regu- 
larization. Here we will present a formal derivation implying the ultraviolet regularization 
(e.g. a lattice regularization). 

In the phase in which electric topological defects condense (while magnetic ones are 




(76) 
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dilute) the partition function requires a formal sum also over the form Qi 



exp 



Z = J VAVBVQ 



(77) 



Let us now compute the expectation value of the 't Hooft operator, {Lh), which represents 
the amplitude for creating and separating a pair of vortices with fluxes ±(p: 



exp 



(Lh) = ^ / VAVBVQ 
If (^1 ^ dBd-i +Al^ *Qi 



B. 



d-l 



Using Stokes' theorem we can rewrite 



(78) 



Bd-i = / dBd-i , (79) 

where the surface is such that dSd = S^-i and represents a compact orientable surface 
on Mj:. Inserting ( 1791 ) in dTSl ) and integrating over the field A we obtain: 



{Lh) oc fVBVQ6idBd^i + *Qi) 



exp 



Integrating over B gives then: 

(Lh) =oc J VQ exp -'^-^4' *Qi 



(80) 



(81) 



The Poisson summation formula implies finally that the 't Hooft loop expectation value 
vanishes for all flux strengths (p different from 



2tt 



n 



n£ N . 



(82) 



This is nothing else than the Meissner effect, illustrating that the electric condensation phase 
is superconducting. Indeed, the electric condensate carries ki fundamental charges of unit 
l//c2 as is evident from dTT] ). and correspondingly vortices must carry an integer multiple 
of the fundamental fluxon 27r/{ki/k2). All other vorticities are confined: in this purely 
topological long-distance theory the confining force is infinite; including the higher order 
kinetic terms (T/T] ) and the UV cutoff one would recover a generalized area law. 

Another way to see this is to compute the current induced by an external elec- 
tromagnetic field Aext- The con^esponding coupling is /^^^^^ Aext A (*ji + *Qi) oc 
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/j^^^^ ^ext A {dBd-i + *Qi)- Since A^^t can be entirely reabsorbed in a redefinition of the 
gauge field Ai, the induced current vanishes identically, jind = 0. This is just the London 
equation in the limit of zero penetration depth. Including the higher-order kinetic terms 
for the gauge fields and the UV cutoff one would again recover the standard form of the 
London equation. 

Associated with the confinement of vortices there is a breakdown of the original U(l) 
matter symmetry under transformations Ai ^ Ai + dX. To see this let us consider the 
effect of such a transformation on the partition function (1771) with an electric condensate. 
Upon integration by parts, the exponential of the action acquires a multiplicative factor 

exp iP- ( [ AA*Qi-/ AA *Qi) . (83) 

Assuming a constant A, we see that the only values for which the partition function remains 
invariant are 

A = 27r n , n = 1 . . . fci , (84) 

ki 

which shows that the global symmetry is broken from U(l) to Z^^. Note that this is not 
the usual Landau mechanism ofspontaneous symmetry breaking. Indeed, there is no local 
order parameter and the order is characterized rather by the expectation value of non-local, 
topological operators. 

The hallmark of topological order is the degeneracy of the ground state on manifolds 
with non-trivial topology as shown by Wen [32]. In (2-1-1) dimensions the degeneracy for 
the mixed Chern-Simons term was proven in [33] for the case of integer coefficient k of the 
Chern-Simons term. 

The degeneracy of the ground state of the BF theory on a manifold with non-trivial 
topology was proven in [9] in (3-1-1) dimensions. This result can be generalized to com- 
pact topological BF models in any number of dimensions [31]. Consider the model (1671) 
with k = ^ on a manifold x Ri , with a compact, path-connected , orientable d- 
dimensional manifold without boundaries. The degeneracy of the ground state is expressed 
in terms of the intersection matrix Mmn [34] with m,n = l....Np and Np the rank of 
the matrix, between p-cycles and (d-p)-cycles. Np corresponds to the number of genera- 
tors of the two homology groups Hp{Md) and Hd~p{Md) and is essentially the number 
of non-trivial cycles on the manifold M^. The degeneracy of the ground state is given by 
\kik2M\^p , where M is the integer- valued determinant of the linking matrix. In our case p 
= (d-1) and the degeneracy reduces to 

\kik2M\^^-^ . (85) 

6 Frustration 

The gauge theory formulation of JJA [12] clearly shows that the superconducting ground 
state is a P- and T- invariant generalization of Laughlin's incompressible quantum fluid. 
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The simplest example of a topological fluid [18] is a ground state described by a low energy 
effective action given solely by the topological Chern-Simons term 

S = k/A-n J (fx A^e'^'^'duAo, 

for a compact U{1) gauge field A^, whose dual field strength = e^^^°^di,Aa yields the 
conserved matter current. In this case the degeneracy of the ground state on a manifold of 
genus g will be (or {kik2Y if /c = ^1/^2 is a rational number). For planar unfrustrated 
JJA one finds that the topological fluid is described by two k = \ Chern-Simons gauge 
fields of opposite chirality and, thus, there is no degeneracy of the ground state [12], [21]. 

In [21] we argued that frustrated JJA may support a topologically ordered ground state 
with non-trivial degeneracy on higher genus surfaces, described by a pertinent supercon- 
ducting quantum fluid, thus providing a more interesting example of a system in which 
superconductivity arises from the topological mechanism proposed in [1] rather than from 
the usual Landau-Ginzburg mechanism. The role of frustration in the determination of the 
ground state degeneracy of JJA (in the limit Ej 3> Ec) was also analyzed in [35]. 

The fact that frustration can change the coefficient of the Chern-Simons term can be 
easily understood by the following example. Let us briefly review the mean field approach 
to spin liquid states [36]. This description is based on the Hubbard model that, follow- 
ing Anderson, is thought to be relevant for high Tc superconductivity. At half filling this 
model reduces to the Heisenberg model. Although very little is known rigourosly about the 
Heisenberg model on a square lattice, it is believed that the Heisenberg antiferromagnet is 
Neel ordered at T = 0. Is it possible to drive the model toward a disordered ground state to 
give rise to a spin liquid? The answer to this question is yes, and the key element to reach 
disorder is frustration. An example is to consider next to nearest neighbor interactions in 
the Heisenberg model: the result is a chiral spin liquid that breaks P- and T-invariance [37]. 
We will thus concentrate on the Heisenberg model, taking the point of view that the var- 
ious spin liquid phases can be realized with small modification of its Hamiltonian [19]. 
Moreover our analysis can be exactly repeted for the Heisenberg model with next to nearest 
neighbour interactions presented in [37]. 

We will study a mean field theory for the Heisenberg antiferromagnet originally pro- 
posed in [36]. To this end we consider a system of spin S = 1/2 on a square lattice with 
nearest neighbor interactions and periodic boundary conditions. The Hamiltonian that de- 
scribes this model is: 

H = Y, J,, SiSj . (86) 

m 

To obtain the mean field ground state for the spin Uquid we use a fermion representation of 
the spin operator: 

Si = -C^ ia(y apCip , (87) 

where a are the Pauli matrices. The operators (|87] ) reproduce the spin 1/2 algebra of only if 
c^iaCia = l^ namely if there is only one fermion per site (half filling). Substituting ( [87] ) in 
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(l86l ) we obtain (up to an additive constant): 

H = - ^ JijC^iaCjaC^ jpCip . (88) 

m 

Introducing the Hubbard-Stratonovich field Xij we can rewrite the four-fermion interaction 
as (we will concentrate here only on this term): 

F = ... - ^ C^iaXijCja ■ (89) 

We will now treat the field Xij within the mean field approximation xlj^ and, since it is 
complex, we will parametrize it with an amplitude pij and a phase Aij. Self consistency 
implies x)f = (c+ioCja). 

Various possible mean-field solutions have been studied. We will be interested in the 
case in which x^j^ generates a flux: the fermions behave as though they would be moving 
in a magnetic field and the P- and T-symmetries are broken; moreover, when an integer 
number of Landau levels is filled, the electron gas is incompressible. For the case in which 
the flux per plaquette is vr, (half of quantum of flux), it has been shown that, after integrating 
out the fermion fields, the effective theory contains a Chern-Simons term with a coefficient 
k = 2 [37](note that in the construction presented in [37] in a certain range of parameter P- 
and T-symmetry are broken even with a flux vr per plaquette). This corresponds to the case 
in which the first Landau level is completely filled. It also been shown that for generalized 
flux phases, where the flux per plaquette is 2T:p/q with q an even integer, the effective 
theory contains a Chern-Simons term with a coefficient k = q (for g = 2 we obtain the 
previous case, namely 1/2 quantum flux per plaquette). In this case we have q/2 Landau 
levels that are filled. 

Let us now consider the hamiltonian: 

H = J2 JijSiUi.Sj , (90) 

(ij) 

where Uij is a frustration field that satisfies : Hp Uij = exp 2i7r/. Here / is the frustration 
parameter. It acts as a fictitious magnetic field, and, due to the invariance of the theory 
under / — >/ + n,nGZ and / —>■ —f we have < / < ^. Moreover since we have 
imposed periodic boundary conditions the flux must be a rational multiple of the quantum 
flux: / = 2iTp/q. Let us follow the same line of reasoning as before. To this end with 
start with a flux vr per plaquette (first Landau level filled), coming from the mean field 
approximation of ([89l). In this case we end up with a flux: 2'k{1/2 + p/ q) = 2TT{q + 2p) /2q 
per plaquette. We will end up with the problem of the hopping of fermions in a fictitius 
magnetic field that is changed from vr 2Tr{q + 2p) /2q, with q Landau levels filled. This 
frustrated model will be described by a Chern-Simons theory with coefficient k = 2q while 
the unfrustrated model corresponds to a Chern-Simons theory with coefficient k = 2. 



Superconductors with Topological Order and their Realization in Josephson Junction 

Arrays 25 



Let us now go back to JJA. It has been suggested that, in presence of rig offset charge 
quanta per site and n<^ external magnetic flux quanta per plaquette in specific ratios, Joseph- 
son junction arrays might have incompressible quantum fluid [2] phases corresponding to 
purely two-dimensional quantum Hall phases for either charges [38] or vortices [14,39]. 

In [12] we have shown that, if quantum Hall phases for charges or vortices are real- 
ized, then JJA naturally support a topologically ordered ground state and a phase in which 
they behave as a topological superconductor [1]; there is, in fact, a renormalization of the 
Chem-Simons coefficient yielding a non-trivial ground state degeneracy on the torus (and in 
general on manifolds with non-trivial topology). To implement a torus topology we impose 
doubly periodic conditions at the boundary of the square lattice. 

To study this degeneracy let us consider the low energy limit of the partition function 
(OTI ). For the following analysis we will take the time as a continuos parameter (we remeber 
that the discretization of time is allowed by the fact that charges are quantized, but is not a 
necessary condition). In this limit we have: 

Z =Y. I^^i^l exp(-5) , 

{Qo} 

S = fdtJ2 A,,K^,B, + iAoQo + iBoMo . (91) 

X 

This form of the partition function holds true also with toroidal boundary conditions. With 
continuos time the operator Kf^i^ is defined by Kqq = 0, K^i = —eijdj, Kiq = Sieijdj 
and Kij = —Sieijdo, in terms of forward (backward) shift and difference operators Si (Si) 
and di (di). Its conjugate K^,y is defined by Kqq = 0, i^oi = —SiCijdj, Kiq = eijdj and 
Kij = —Sjeijdo. 

The topological excitations are described by the integer-valued fields Qo and Mq and 
represent unit charges and vortices rendering the gauge field components Aq and Bq in- 
tegers via the Poisson summation formula; their fluctuations determine the phase dia- 
gram [12]. In the classical limit the magnetic excitations are dilute and the charge exci- 
tations condense rendering the system a superconductor: vortex confinement amounts here 
to the Meissner effect. In the quantum limit, the magnetic excitations condense while the 
charged ones become dilute: the system exhibits insulating behavior due to vortex super- 
conductivity accompanied by a charge Meissner effect. 

By rewriting the topological excitations as the curl of an integer-valued field 

Qo =Ko^Yi, YieZ, 

Mo =KoiXi, Xi^Z, (92) 

we get the mixed Chem-Simons term as follows: 

Z = E{x.,yo fVA^fVB^ exp (S) , 
S =-i^iIdtE^ AoKoi {B, - 2T,Yi) 

+Bokoi {Ai - InXi) + AiK.jBj . (93) 
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^From (1931 ) one sees that the gauge field components Ai and Bi are angular variables due to 
their invariance under time-independent integer shifts. Such shifts do not affect the last term 
in the action, which contains a time derivative, and may be reabsorbed in the topological 
excitations Xi and Yi, leaving also the first term of the action invariant. The low energy 
theory is thus compact. 

In analogy with the conventional quantum Hall setting one should expect the charge 
and vortex transport properties to depend on the ratios of the offset charges (i.e. the filling 
fractions) {nq/n^j,) and {n^/nq), respectively. Due to the periodicity of the charge- vortex 
coupling, however, (n^) is defined only modulo an integer as far as charge (vortex) 
transport properties are concerned. Using this freedom one may define effective filling 
fractions (we shall assume > 0, n,^ > for simplicity) as 

Vr, = "1 ■ , < < 1 , 



[n4,\ +[nq\ • ^ 
~i-<l> 



^ r r 1+ ' 0<z.<^<l, (94) 

where indicate the smallest (greatest) integer greater (smaller) than nq. Of course, 

these effective filling fractions are always smaller than 1. 

In [12] we assumed the existence of these quantum Hall phases and discussed them in 
the framework of the gauge theory representation of Josephson junction arrays, showing 
that, depending on certain parameters of the array there are both a charge quantum Hall 
phase and a vortex quantum Hall phase. Here we will concentrate on the low energy limit 
of the charge quantum Hall phase and we will show that the system has topological order 
and behaves as a superconductor when charge condenses. 

The pertinent low energy theory is now given by: 

S= [dtV] --A^K^^B, - 'J^A^K^^A, , (95) 
J ^ IT vr 

with Uq = p/n. The main difference with ( |9T] ) is the addition of a pure Chern-Simons term 
for the A^ gauge field . We have also rescaled the coefficient of the mixed Chem-Simons 
coupling by a factor of 2 (compare with (|9T| ) ). This factor of 2 is a well-known aspect of 
Chern-Simons gauge theories [40]. Moreover, since in JJA the charge degrees of freedom 
are bosons, the allowed [12] filling fractions are given by f ^ = ^ , with pn = even integer 
in accordance with [41]. As a result, the action ( |95] ) may now be written in terms of two 
independent gauge fields A^ and B'^ = B^ + I'qA^ yielding: 



S 



J dtJ2-^\K^.B^,. (96) 



In describing JJA one has to require the periodicity of charge-vortex couplings; the 
coupling of the topological excitations enforcing the periodicity of the mixed Chern-Simons 



term A^K^j^i^Bf, is then: 



S = J dtJ2--- + ipAoQo + inBoMo , (97) 
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that can be rewritten as: 

S = J2--- + i^pAo (Qo + Mo) + UnB^^Mo . (98) 

X 

Due to the replacement ^ B'^, the periodicities of the two original gauge fields are 

Ai Ai + Tin ai , ai £ Z , 

Bi ^Bi + npbi, keZ, (99) 

and 

Bl -^Bf + TTpbi, keZ . (100) 

The resulting low energy theory is thus, again, compact. 

Using the representation ( |92l ). one may rewrite the mixed Chem-Simons term as 

X 

_,(eM^Ko.(M._2vrX.) . (101) 

In this representation it is clear that the topological excitations render the charge-vortex 
coupling periodic under the shifts 

2Ai 

A, = ^ A- + 27r , tti e Z , 

n 

b[ = ^ -^B[ + 2^hi, heZ . (102) 

P 

This model corresponds to two Chem-Simons terms with coefficients ^k/Air with k = 
np/2 an integer. It is worth to point out that, since B^ does not have a definite parity (is a 
liner combination of a vector and a pseudovector) the model is not P- and T-invariant, as it 
must be due to the presence of the Chem-Simons term for the field A^. 

The hallmark of topological order is the degeneracy of the ground state on manifolds 
with non-trivial topology [3]. The torus degeneracy on the lattice of the Chern-Simons 
model was computed in [42]. For a single Chern-Simons term this degeneracy is {k)^ where 
k is the integer coefficient of the Chern-Simons term, and g the genus of the surface. In our 
case this degeneracy is 2 x {ky = 2 x since we have two Chem-Simons terms. This 
degeneracy is exactly what is expected for a doubled Chem-Simons model [10], for which 
the physical Hilbert space is the direct product of the two Hilbert spaces of the component 
models. 

We will now demonstrate that the phase where topological excitations Qo condense 
while Mq are dilute describes an effective gauge theory of a superconducting state. The 
partition function is: 



Zle = E / / exp i-S) , 
{Qo} 

/■I ■ 1 

dt A^K^^Bl + ^Ao(27rQo) . (103) 
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To this end note first that a unit external charge, represented by an additional term 
i2-Kao{t, x)5xxo is completely screened by the charge condensate, since it can be absorbed 
into a redefinition of the topological excitations Qq. In order to characterize the supercon- 
ducting phase we introduce the typical order parameter namely the 't Hooft loop of length 
T in the time direction: 

Lh = exp J ^ , (104) 

where (/>o(t, x) = {6 {t + T/2) - 9 {t - T/2))5^^^ - {0 {t + T/2) - 9 {t - T/2)) S^^^ and 
r/2, x), 0j(r/2,x) are unit links joining xi to X2 and X2 to xi at fixed time and 
vanishing everywhere else. Its vacuum expectation value {Lh) yields the amplitude for 
creating a separated vortex-antivortex pair of flux cj), which propagates for a time T and is 
then annihilated in the vacuum. 

Since we replaced B^, we may rewrite the 't Hooft loop as: 



Lh = exp (^ict>^ j dtY, {^f^Bf, - ^^/.'A^)) • (105) 



To compute (Lh) one should integrate first over the gauge field B'^ to get 

{Lh) oc ^ / VAf, 6 (k^^A^ - ^cp^) x 
{Qo} 

expic/.^ [^j dtY, ^o(2^Qo) - ^M^j ■ (106) 



The sum over Qq enforces the condition that -^Aq be an integer. As a consequence, 
KiqAq = = (j)(f>i with rij an integer. We thus have: 

cP=^q,qGZ- (107) 

thus, {Lh) vanishes for all fluxes different from an integer multiple of the fundamental 
fluxon, which is just the Meissner effect. In the low-energy effective gauge theory vortex- 
antivortex pairs are confined by an infinite force which becomes logarithmic upon including 
also higher-order Maxwell terms. 

^From (IIOII ) one has (depending if p is an even integer or n is an even integer) either 
that n is the charge unit with p/2 units of charge, or, viceversa, that p is the charge unit 
with n/2 units of charge. By rewriting (11071) as = ^27rn one finds the standard flux 
quantization (j) = j^n where e is the charge unit and N is the number of units of charge. 

In this paper we have derived a superconductivity mechanism which is not based on the 
usual Landau theory of spontaneous symmetry breaking. Our considerations here focused 
on the low-energy effective theory in order to expose the physical basis of the topologi- 
cal superconducivity mechanism. It is however crucial to stress that the simplest example 
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(k=l) of this type of topological superconductivity is concretely realized as the global su- 
perconductivity mechanism in planar Josephson junction arrays, as we have shown in [12]. 
Naturally, it would be most interesting to find examples of microscopic models realizing 
this superconductivity mechanism with more complex degeneracy patterns. A possibilty is 
represented by models with frustration. The frustrated planar JJA we have analized provide, 
in the quantum Hall phases, an explicit example of both topological order with non-trivial 
ground state degeneracy on manifolds with non-trivial topology, and of a new supercon- 
ducting behavior [21] analogous to Laughlin's quantum Hall fluids. 
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